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Abstract 

Based on the previous work in [Y. Jiang, Y.M. Shi, H.T. Feng, W.M. Sun and H.S. Zong, Phys. 
Rev. C 78, 025214 (2008)] on the quark-meson vertex and pion properties at hnite quark chemical 
potential, we provide an analytical analysis of the weak decay constant of the pion (/tt[^]) and the 
pion mass (m n [u}) at finite quark chemical potential using the model quark propagator proposed in 
[R. Alkofer, W. Detmold, C.S. Fischer and P. Maris, Phys. Rev. D 70, 014014 (2004)]. It is found 
that when /x is below a threshold value /xo (which equals 0.350 GeV, 0.377 GeV and 0.341 GeV, for 
the 2CC, 1R1CC and 3R parametrizations of the model quark propagator, respectively.), f n [fj] and 
m-vd/i] are kept unchanged from their vacuum values. The value of /io is intimately connected with 
the pole distribution of the model quark propagator and is found to coincide with the threshold 
value below which the quark-number density vanishes identically. Numerical calculations show that 
when u becomes larger than uq, /^[/i] exhibits a sharp decrease whereas m-^lfj] exhibits a sharp 
increase. A comparison is given between the results obtained in this paper and those obtained in 
previous literatures. 
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The in-medium modification of the properties of the pion is of fundamental interest in 
hadron physics. The pion is identified as a Goldstone boson arising from the spontaneous 
breakdown of chiral symmetry which is essential for describing low-energy hadronic phenom- 
ena. Since chiral symmetry is expected to be restored at high enough density, the change 
of pion properties in medium will provide crucial information on the restoration of chiral 
symmetry. Among these, the weak decay constant of the pion f n and the pion mass are 
the two most important quantities, since they are closely related to the spontaneous break- 
down of chiral symmetry of Quantum Chromodynamics (QCD). Unfortunately, so far it has 
not been possibile to obtain detailed information about modification of pion properties in 



medium directly 
sort of problems 



rom QCD. In this situation, different models have been used to study this 



111 ], the pion 



J 13, S, fl E, 3, 3, i Q. Just as was pointed out in Ref. 
has a dual role: it can be identified as a quark-antiquark bound state as well as a Goldstone 
boson arising from the spontaneous breakdown of chiral symmetry. From the point of view 



that the pion can be regarded as a quark-antiquark bound state, the full dynamical informa- 
tion of the pion is contained in the corresponding Bethe-Salpeter Amplitude (BSA): T w (k,p) 
(k is the relative and p the total momentum of the quark-antiquark pair), which is the one- 
particle-irreducible, fully-amputated quark-meson vertex. The Dyson-Schwinger equations 
(DSEs) of QCD provide a nonperturbative, continuum framework for analyzing such quark- 



meson vertices directly 
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fl[ 



, ll5( . The aim of this paper is to study the change of 
fn and m n with quark chemical potential ii in the framework of this nonperturbative QCD 
model. 

The DSEs of QCD have been used extensively at zero temperature and zero quark chem- 
ical potential to extract hadronic observables [l2,[l^, Q, 3|- However, this is very difficult 
at finite quark chemical potential due to the fact that the number of independent Lorentz 
structures of the quark-meson vertex at finite /i is much larger than that of the corresponding 
one at \x = 0. In Ref. 16|, using the method of studying the dressed quark propagator at 
finite jx given in Ref. [13] , the authors have given a new approach for tackling this problem. 
Based on the rainbow-ladder approximation of the DSEs and the assumption of analyticity 
of the quark-meson vertex in the neighborhood of 11 = and neglecting the /i-dependence 
of the dressed gluon propagator, the authors show that the general quark-meson vertex 
at finite ii can be obtained from the corresponding one at \x = by a shift of variable: 
T\p](k,p) = T(k,p), where k = {k^k^ + ip). From this result the authors of Ref. lf| 



2 



numerically calculated fn\p>] and m^-[/x] for // < 300 MeV. It is found that f w \p] increases 
slowly (with an increase of less than about 0.01%) and m^-[/x] falls slowly (with a decrease 
of less than about 0.06%) with increasing //. Numerically the change of f w \p] and m^/i] 
is so small that one can think f n \p] and m T [/i] does not change with \x for /i < 300 MeV 
within numerical errors. One of our motivations for this work is to explore the mathematical 
reason behind this. Based on the work in 16}], in this paper we provide an analytic analysis 



of ./*■[/*] and m^f/i]. It is found that when p, is below a critical value //o, fn\p] and m^^] 
are kept unchanged from their vacuum values. Moreover, numerical calculations show that 
when p becomes larger than p , f n [tA exhibits a sharp decrease whereas m^f/i] exhibits a 
sharp increase. 

According to Ref. }16| . the pion decay constant at finite p can be expressed as the 
following 



tr 



T 



-jl5luS(q + )Ti(q;p)S(q. 



(1) 



where S(q) is the full dressed quark propagator, q± = q ± p/2, q = {q.q^ + ip), y are 
the flavor SU(2) generators and / = / d 4 q/ \2ir) 4 . In the present paper we will not write 
the renormalisation constants explicitly because one would find that in the final result the 



renormalisation constants cancel each other. In fact, Eq. (JTJ) is the expression o 



f n which 



is independent of the renormalisation point and the regularisation mass-scale 11]. 
The integral of the right-hand-side of Eq. (CQ) can be rewritten as: 

d A q _ f d 4 q 



(27T) 



(27T) 



(2) 



Contracting both sides of Eq. (TjQ) with p u and using Eq. (J2J), we obtain the following: 



5 ij fM 
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oo+ifi 
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p 2 J (2vr) 3 J ox (2tt) 



■tr 



■tr 



T 



-ls^S(q + )Ti(q;p)S(q-) 



lB^S(q + )Tl(q;p)S(q- 



(3) 



where the integration path C\ is depicted in Fig. 1. 
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FIG.l. The integration path in the complex plane. 
Let us use z n = Xn + i^n (w n > 0), n — 1, 2 • • • to denote the poles of the function 



1 



—^S{q+)T 3 n {q;p)S{q. 



(4) 



located in the upper half complex g 4 plane. According to Cauchy's theorem we obtain the 
following from Eq. ([3]): 

d 3 q 



dq A 



p 2 J_ (2tt) 3 Jcj (2tt) 
1 f d 3 q f dq A 



p 2 J_ (2tt) 3 Jc (2tt) 
d 3 q 



tr 



tr 



r 



— 7 5 /S'(g + )r^(9;p)5'(g- 



-oo 



T 



-^ n f5$S(q + )Tl(q;p)S(q- 



J j^kY, d (»-"n)Res{FV(zy,z n } 



(27T) 3 ^ 



(2vr) 3 ^ 



^^(/x-^ n )Res{F^(z);^}. 



(5) 



From Eq. (jHJ) it is easily seen that when fi < min{uj n }, the function F^^) has no pole 
in the region Q (the region enclosed by G\ and C , see Fig. 1) and therefore f n [fj] = f n , 
which means that for small enough \x the pion decay constant should be independent of \x. 
Of course, when \i > min{uj n } the pion decay constant can have an explicit /i-dependence. 

In the chiral limit, expanding the trace term of the right-hand-side of Eq. (j4]) to 0(p 2 ) 
near p = [3] , we have the following: 



F ij {q 4 



i r t 1 



P 



2 -i 2 



5 + -p • 95 



ri(g,o) + o( P ) 7£ 



5 - -p • 95 



(6) 



where we have adopted the approximation 12] 



Ti(q,p) = Ti(q,0) + O(p) l5 . 
With this approximation T 3 w (q,0) can be expressed as [12I. fis| 



ri(g,0) = r^5 



(7) 



(8) 



where B(q 2 ) is the scalar part of S 1 (q). Noticing that ti^p'S-y^S] = 0, we obtain the 
following: 



p-dSTi(q,0)S-STi(q,0)p-dS 



+ 0{p). 



(9) 



Substituting Eq. (jSj) into Eq. (jUJ) and using tr(rV- J ) = 25 u , we obtain 

p • d-S^S - S75P • dS 



Adopting the following expression of S(q) 

1 



5(g) 



i^4((r«) + 5(g») 



-i{ja v (q 2 ) + a s (q 2 ), 



we obtain 



F ij {q, 



1 h,- 1 8(7, 



2 <t„V + a 2 



, 2(p-g)% , , 
c s Cu H 5 — (crs^ - cr s er„J 



In 



where ' means d/dq 2 and 



4a, 



Then Eq. (j3j) can be written as 



(p-qf 



a s a v + 2- 



p* 



(10) 

(11) 

(12) 
(13) 

(14) 



U\p] ^f K - — J 77^31 E 9 (V ~ ^n)Res{F(z); z n }. 



f* L (2vr) 3 ^ 



(15) 



To determine the pole distribution of function F{gA, we should first specify the form of 
the dressed quark propagator. Here, as in Refs. |16|, ll9] we adopt the following propagator 
proposed in Ref. [201 ]: 

np / 



S(q) = E 



T 3 



+ 



(16) 
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The propagator of this form has np pairs of complex conjugate poles located at aj ± ibj. 
When some bj is set to zero, the pair of complex conjugate poles degenerates to a real pole. 
The restrictions of the parameters rj, aj and bj in the chiral limit are 20| 

up 



E r i = 

3=1 

Hp 

E T 3 a 3 



1 
2' 



0. 



(17) 
(18) 



If we are not in the chiral limit, the right hand side of Eq. (fT8"|) should be replaced by the 
current quark mass. The value of these parameters are shown in Table I, where 2CC, 1R1CC 
and 3R stand for three meromorphic forms of the quark propagator, respectively: two pairs 
of complex conjugate poles, one real pole and one pair of complex conjugate poles, three 
real poles. 

Table I. The parameters used in the calculation of F(qi) and f n . These parameters are taken directly from Ref. [20I ] . 



Parameterization 


n 


ai (GeV) 


61 (GeV) 


T2 


a 2 (GeV) 


62 (GeV) 


r-3 


a 3 (GeV) 


2CC 


0.360 


0.351 


0.08 


0.140 


-0.899 


0.463 






1R1CC 


0.354 


0.377 




0.146 


-0.91 


0.45 






3R 


0.365 


0.341 




1.2 


-1.31 




-1.06 


-1.40 



Without losing generality we assume p v = (0,p) (i.e. the pion is at rest) and write 



(p-q) 2 q\p 2 



(19) 



p* p± 

Now let us calculate F{q^). With the quark propagator given in Eq. (1161) we can obtain 



F{qA 



U[q 2 + (a, + ibj) 2 ] 2 [q 2 + (% - ibj) 2 } 2 U(q 2 + vl) ' 



(20) 



3 k 

where 5 is a polynomial of q\ (for the detailed calculation of F(q4), H and rjk, see the 
Appendix). The values of % are shown in Table II (% are ordered from small to large 
according to their real part). 

Table II. The calculated values of r\^. 



Parameterization 


m (GeV) 


m (GeV) 


m (GeV) 


m (GeV) 


V5 (GeV) 


2CC 


0.350 


0.723-0.351i 


0.723+0.351i 






1R1CC 


0.377 


0.723-0.328i 


0.723+0.328i 






3R 


0.341 


0.617 


1.31 


1.40 


1.849 



Here it should be noticed that when some bj = (the quark propagator has a real pole), 
some r/k must exactly equal the corresponding \a,j\ (see the Appendix). For 1R1CC case, 



b\ = and 771 = \ai\. For 3R case, all bj = and 771 = |ai|, 773 = | a-2 1 , Va = | a-3 1 - For 2CC 
case, because 61 = 0.08 GeV is very close to zero, the value of 771 is very close to a\. 

Because q 2 = ql + q 2 , according to Eq. (I2"0"j) the poles of Fi^q^): z n = Xn + iu n are decided 
by the following equation 

( X n + lUn? + ? + {inR + = 0, (21) 

where and £ n j are the real and imaginary part of r/k or a,j ± ibj . One can easily find 



n ~ \ 2 [ ' 

Xn = • (23) 

From Eq. (I2"2"j) we find that for /ii < \^ u r\ the corresponding uo n is always larger than fi, 
irrespective of q. For /i > \£ nR \, u n < /i when q 2 < /i 2 - (^^///i 2 ) - £ 2 R + £ 2 /, and u n > 
when (f > fi 2 — (£ 2 /j£ 2 / / '/i 2 ) — £ 2 R + £ 2 /. Therefore Eq. ( TloT) can be written as 

A„(/i) 

AM = ^- 9 -VE^-|e^l) / ^Res^);^}, (24) 

^ Jit n J 

where 



An(/i) = vV - (^e 2 //^ 2 ) - a R + e n i ■ (25) 

From Eq. ( 1241 and the values of a.,-, 6j and 77^ in Table I and II we find that when [i is below 
some threshold value fio, the pion decay constant at finite chemical potential f w \p] is kept 
unchanged from its vacuum value. The threshold value /zo, which equals the minimum of 
the real part of <%,• ± ibj and 77^, is shown in Table III. 

Table III. The calculated values of /iq. 



Parameterization 


Ato (GeV) 


2CC 


0.350 


mice 


0.377 


3R 


0.341 



Here we note that in Ref. [2jJ it is found that when /i is below the same threshold value 
fiQ, the quark- number density vanishes identically. Namely, fi = fiQ is a sin gula rity which 



separates two regions with different quark-number densities. In fact, in Ref. 22|, based on 



a universal argument, it is pointed out that the existence of some singularity at the point 



H = Ho and T = is a robust and model-independent prediction. Below // = /iq, the QCD 
system at finite /i remains in the vacuum (ground state) of QCD at /i = 0, so the properties 
of the Goldstone boson excited from this vacuum does not change with /i. Thus the result 
that f n \pi] is kept unchanged from its vacuum value is just to be expected. Here it should 
also be noticed that in our method the value of /io is intimately connected with the pole 
distribution of the quark propagator. 

In Ref . [3] , with the same quark propagator the authors find that the quark condensate 
at finite chemical potential is kept unchanged from its vacuum value when /i < /i . From the 
Gell-Mann-Oakes-Renner relation fl[n]m\[ji\ = 2m(qq) [fi] + O(m 2 ) [111, ll6j (where m^/i] is 
the pion mass at finite /i, m is the current quark mass and (qq)o[fj] is the quark condensate 
in the chiral limit at finite /i) one would also conclude that m n \p\ is kept unchanged from its 
value at = when /i < /io- In Ref. [la ], the authors did not made an analytical analysis 
of f n [fj] and m,^] by the method of pole analysis, but instead made a direct numerical 
calculation. There exist numerical errors in this calculation. Within numerical errors /7r[/i] 
and m,r[/x] do not change with \i for /i < 300 MeV. The analytical analysis made in this 
paper explains the numerical results obtained in 16]. 

For jj, > fi , one can calculate f n \p] and m^f/i] numerically based on Eq. (jUj) and the 
Gell-Mann-Oakes-Renner relation. The behaviors of f n [fj] and m^/i] for // > /io are shown 
in Fig. 2 and Fig. 3. One sees that f n [fi] exhibits a sharp decrease whereas m„-[/i] exhibits 
a sharp increase near /io for all three cases. This result is quite different from the result 
in previous literatures. For example, in a recent work those authors also investigated 
f n and m n at finite density within the framework of the nonlocal quark model from the 
instanton vacuum. Their results show that in the range < /i < 320 MeV, f v falls slowly 
whereas m n increases slowly. This behavior of and m n is qualitatively different from that 
found in this paper. 

Finally, we should emphasize that in obtaining our results about f n [fj], m n [fj] and {qq)o[fj] 
in this paper, we have made these approximations and assumptions: (1) we adopt the 
rainbow-ladder approximation of the DSEs; (2) we assume the quark propagator and quark- 
meson vertex are analytic in the neighborhood of /x = 0; (3) we have neglected the /i- 
dependence of the dressed gluon propagator, (for a discussion about these approximations 
and assumptions, see Ref. 16|]). For further study one should consider improvements on 
these approximations. 
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FIG. 2. The \i dependence of /„- near /i . 
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FIG. 3. The \i dependence of near fi . 



To summarize, based on the previous work in Ref. 16| on the quark-meson vertex and 
pion properties at finite quark chemical potential, we provide an analytical analysis of the 
weak decay constant of the pion (/^[/i]) and the pion mass (m n \ u] ) at finite quark chemical 
potential using the model quark propagator proposed in Ref. 20|. It is found that when 
/i is below a threshold value fi (which equals 0.350 GeV, 0.377 GeV and 0.341 GeV, for 
the 2CC, 1R1CC and 3R parametrizations of the model quark propagator, respectively), 
f n [fi] and [/x] are kept unchanged from their vacuum values. The value of /io is intimately 
connected with the pole distribution of the model quark propagator and is found to coin- 
cide with the threshold value below which the quark-number density vanishes identically. 
Numerical calculations show that when /i becomes larger than /io, fir[t~i] exhibits a sharp de- 
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crease whereas m n \p] exhibits a sharp increase. These results are quite different from those 
obtained in previous literatures. For example, our results are qualitatively different from 
those reported in Ref. 7|, which uses the nonlocal chiral quark model from the instanton 
vacuum to investigate f n and m n at finite density. 
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APPENDIX A: THE ANALYSIS OF THE POLES 
1. General Analysis 



With the quark propagator given by Eq. (TIB"]) one can find the following 



E 

3 

E 



+ 



q 2 + (dj + ibj) 2 q 2 + (a, - ibj) 2 j 



k 
h 

U 

q 2 + (dj + ibj) 2 ' q 2 + (a,- - ibj) 2 \ f 



rj{aj + ibj) r ! (a J ib,) 



with 



A = E 2r i (V + S 2 " b d IW + (afc + rt k ) 2 } [q 2 + (a k - ib k ) 2 } 
f s = Y. 2r i a j + a ) + & ?) II il 2 + ( a k + *h) 2 } [q 2 + (a k - ib k ) 2 } 

3 k^3 

j 



Then one has 



F{q, 



4(7., 



a s a v + 2ql(a s a' v — a' s a v 



1 



fof 2 q 2 + f 2 ' 



(Al) 
(A2) 

(A3) 
(A4) 
(A5) 

(A6) 



where 

S = 4f s [fsfv + 2(q 2 - f)(f.f v - f'j v )] . (A7) 

For convenience let us use x 2 = q 2 with x a complex number. Then the denominator of the 
right-hand-side of Eq. (1A6j) can be decomposed as 



fy + fs = (fvX + if s )(f v X-if s ). 



(A8) 
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/„ and f s can be expressed as 



fv = E' 

3 

fs = £■ 



/o 



+ 



fo 



X, 



+ (oj + ibj) 2 x 2 + (aj - ibj)' 



foicij + ibj 



+ 



fo(aj - ibj 



x 2 + (a,j + ibj) 2 x 2 + (cij — ibj) 2 



so one has the following 



(fvX + ifs)(fvX-ifs) 

sr^ <• \ x + i(dj + ib~) x + i(a,j — ibj) 
E ''/•/(' > ■>■ - 

x |Eo/o 

E r j/o 



x 2 + (dj + ibj) 2 x 2 + (dj — ibj) 2 
x — i(dj + ibj) x — i(dj — ibj) 



x 2 + (dj + ibj) 2 x 2 + (dj — ibj) 2 



1 



+ 



1 



x < E '•//o 



x — i(dj + ibj) x — %{dj — ibj) 
1 1 



+ 



x + i(dj + ibj) x + i(dj — ibj) 



fo can be expressed as 



fo = T[[ x + i ( a ki +ib kl )][x + i(d kl -ib kl )] 
fci 

x Y[[x-i(d kl + ib kl )][x - i(d kl -ib kl )]. 
k 2 



Therefore one obtains 
E r J'/o 



+ 



(A9) 
(A10) 



(All) 



(A12) 



x — i(dj + ibj) x — i(dj — ibj) 
= E \j'A x ~ ^( a i ~ ibj) + x - i(dj + ibj)] Y[[x + i(a kl + ib kl )][x + i(a kl - ib kl )\ 



3 " ki 

x Yl [x - i(a k2 + ib k2 )}[x - i(d k2 - ib k2 )] 

k2^j 

Y[[x + i(d kl + ib kl )][x + i(d kl - ib kl )] 



x E 2 ^^ - ia j) ]l[ x - K a k + ib k )][x - i(d k - ib k )\ 

3 k^j 



(A13) 



and 



E r i/o 

3 



+ 



x + i(dj + ibj) x + i(dj — ibj) 
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= £ \ r j \ x + *( a i - %) + x + ^( a i + II l x + ^( a ^i + *fci)H x + ^( a fci _ *fcj] 

x J|[x - i(a k2 + ib k , 2 )}[x - i(a k , 2 - ib k2 )} 
= Y[[x - i(a k2 + ib k2 )][x - i(a k2 - ib k2 )) 

xJ2 2r A x + ia j) T[l x + i{a k + ib k )}[x + i(a k - ib k )\. (A14) 

3 k+j 

With Eq. flA12j) one can find the following 



fv X + fs = (fvX + if s )(f v x -if a ) 

= fo I 5Z 2r i( x ~ ia i) Iii x - *( a fc + ib k)}[x - i{a k - ib k )} 
I j k+i 

x \ 5Z 2r i( x + ia j) T[[ x + ^ a k + ib k )\[x + i(a k - ib k )} > . (A15) 

{ j k+j J 

Hence, in order to determine the poles of F^q^), one should solve the following three equa- 
tions: 

fo = II^ 2 + ( a i + ib j)V + (a, - tbj) 2 } = 0, (A16) 
j 

J22rj(x -idj)Y[[x -i(a k + ib k )][x -i(a k -ib k )] = 0, (A17) 

j k^j 

^2r j (x + icij) Y[[x + i(a k + ib k )][x + i(a k - ib k )] = 0. (A18) 

j k^j 

Here it should be noted that if some bj = then x = iaj (or x = —idj) must be the solution 
of Eq. (1A17I) (or Eq. (1A18P ). One should also be aware that after finding the roots of the 
above equations one should substitute them into H to ensure that S(x) ^ (we will see it 
in the discussion of 1R1CC and 3R case below). For general rip Eq. (1A17|) (or Eq. (1A18H ) 
is an equation of degree 2np — 1 in i and it is almost impossible to give the analytic form 
of the solution for general rj, cij and bj when np > 2. 

2. Detailed Calculation of the Poles 

For 2CC case one can find the following 

f v = q 6 + d vl q 4 + d v2 q 2 + d v3 (A19) 
/* = d sl q 4 + d s2 q 2 + d s3 (A20) 
fo = [q 4 + 2(a? - bl)q 2 + (a\ + bj) 2 ]^ + 2{a\ - b\)q 2 + (o| + b\)\ (A21) 
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where d v i, d v2 , d v3 , d s i, d S 2, d s s are coefficients decided by rj, a,j, bj. With parameters shown 
in Table I the solutions of Eq. (IA17I) and Eq. (IA18I) are found to be rjx = 0.350 GeV, 772,3 = 
(0.723 ± 0.351?) GeV. Of course, one can directly verify 

fW + f s = k{q 2 + rfttf + VIW + Vl). (A22) 
So the poles of F(q^) for 2CC parameters (in the upper half complex g 4 plane) are 



Zl 


— i\J Q 1 + Vi (simple pole) 


(A23) 




= X2 + iu> 2 (simple pole) 


(A24) 


Z3 


= X3 + ^3 (simple pole) 


(A25) 




= Xi + ^4 (double pole) 


(A26) 


Z5 


= Xs + (double pole) 


(A27) 


Z6 


= xe + iwe (double pole) 


(A28) 


Z7 


= Xi + ^7 (double pole) 


(A29) 



with 



X2 



\ 



q 2 + (Rer] 2 ) 2 - (lmr] 2 ) 2 + J [q 2 + (Rer] 2 ) 2 - (lmr] 2 ) 2 ] 2 + A(Rer] 2 ) 2 (lmr] 2 ) 2 



-X3 



(Rer/ 2 )(Imr/ 2 ) 

<jj 2 



X4 = -X5 

UJq = IO7 = 
X6 = ~X7 



? + 


-b 2 + 


+ a? 


- b 2 ) 2 + Aa\b\ 


i 




2 




a\bi 








004 








q 2 + a 2 2 


-63 + 


+ 4 


- b 2 2 ) 2 + Aa 2 b 2 



\ 



a 2 b 2 
u 6 



(A30) 
(A31) 

(A32) 
(A33) 

(A34) 
(A35) 



For 1R1CC (and 3R) case, the analysis is similar except a little modification for correctly 
analyzing the degree of the poles. Because 61 = for 1R1CC case (for 3R case, all bj equal 
zero) the function f v and f s has a factor of q 2 + a 2 (see Eq. ( 1A3I) and Eq. (I A4I) ) which 
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would be canceled by the same factor in f . Therefore for 1R1CC case one should adopt the 
following modified expressions 



fvl 



fs 



q 2 + a\ 

2n[q 2 + (a 2 + ib 2 ) 2 } [q 2 + (a 2 - ib 2 ) 2 } + r 2 (q 2 + a 2 ) [q 2 + (a 2 - ib 2 ) 2 } 

+r 2 (q 2 + a\) [q 2 + (a 2 + ib 2 ) 2 } (A36) 

fs 

q 2 + a\ 



2r x a x \q 2 + (a 2 + ib 2 ) 2 ][q 2 + (a 2 - ib 2 ) 2 ] + r 2 (a 2 + ib 2 )(q 2 + a\)[q 2 + (a 2 - ib V2i 



■2 J 



fi 



+r 2 (a 2 - ib 2 )(q 2 + a\)[q 2 + (a 2 + * 2 ) 2 ] (A37) 
/o 



g 2 + af 



= (g 2 + a 2 )[? 4 + 2(a 2 - 6 2 )g 2 + (a 2 + 6 2 ) 2 ]. (A38) 
According to the decomposition in Eq. ( 1A15I) one has 

flq 2 + f s = fo(q 2 + Vl)(q 2 + %W + vl) (A39) 
with r/i.2,3 being obtained by solving Eqs. ( 1A17I) -( IA18I) . Then 

pfn N = ±Wsfv + 1{?-?)U.f' V -fvf.)\ , A4m 

W4j /o 2 (g 2 + ^i 2 )(g 2 + % 2 )(g 2 + % 2 ) 1 J 



+ 2(g 2 - gK/sLgl - /.l/ jx)] 

/! 2 (g 2 + ^)(9 2 + % 2 )(? 2 + % 2 ) 



(q 2 + ai). (A41) 

Because 771 equal a\ exactly, one would find that iJ q 2 + a 2 is a double pole. The analysis 
for 3R case is similar. 
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